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METHOD FOR BISPECTRAL: PTr K ING OF ANELL IP T I CAL NMO / 
- CORRECTION PARAMETERS 3 - . .. 

The field of this invention is seismic prospecting 
The invention relates particularly. . to. processing of 
seismic traces in a common midpoint gather. 

.., . The . invention , is more precisely related to a method 
of determining velocity y, and anellipticity n parameters 
: ^cessary to perform ^processi^ - correction 

, pf, seismic traces.-. . - . ~: ztr-^. • ~ - * 1 - • - ' .. 
Seismic prospecting usually • consists , of emitting' 
seismic waves in the' subsoil using one or several seismic ' 
sources, and making^ : surface. ; records « of V" seismic * -dat a 
corresponding to seismic waves , reflected . on geological 
; interfaces^, in the f subsoili, ^ 

..pa^ed, reflectors) . using: receivers, (alsoT called geophones 
.or hydrophones), depending on .whether thei prospecting , is 
bexng done on land or, at sea) and then- processing thesW 
data, to extract useful information about the geology of 
the subsoil. 

The record of the seismic energy produced by each 
receiver during data acquisition is called a seismic 
20 trace. 
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One conventional technique for seismic prospecting 
is multiple coverage, in which sources and receivers are 
arranged such that several seismic traces are grouped 
together at the same midpoint (in other words the point 
5 at equal distance between the source and the receiver for 

a given trace) . . ... .-' ■ ■ 

Although' seismic traces contain useful information 
about seismic reflections and" geology of the subsoil, 
they also contain noise components. 
10 One- of the^ firsts purposes! of the . processing of 

seismic data is > tb; v 'eiimihate}:. ^ at . least attenuate., 
these unwanted noise components such that the useful 
information catf be clearly identified and interpreted. 

One conventional method used to attenuate these 
15 noise components is Common MidPoint (CMP) gather . Traces 
with the same midpoint are then grouped as a function of 
, the i: distance between the source and the receiver (called 

vthe'offset) ^ ' •• A-f$'-^ ; Z '^—■{''f^ .' -j/.;-'-" : "' '•"■•.; • 

•-' * In general / representation of seismic data- in images 

20 requires the use of, processing including:" * ^ ; * 

; ' _ a so-called.v'/Tzdv (Transform"^ to" Zero -Offset) , 

: operation' designed. ^.compensate for- the NMO effect' 'of 
. ; paths); by ] correcting 

.the' arrival t imeV of ^traces ; wi thV zero offset, ■ , '••• : -. ; - \ 
25 , ■ - and a migration- operation designed to reproduce 
■•• •• the correct shapes of the geological interfaces. ■; K- 

. . Although: these TZO and:: 'migration- operations ?, are 
usually carried out one after the- other, they can- be 
carried out jointly. This is the case particularly when a 
time migration is done before stack (PSTM - Pre-Stack 
Time Migration) . 



30 



WO 2005/008291 >v 



PCT/TB2004/002617' 



in a simplified manner, the TZO operation simulates 
acquisition of seismic data by sources and receivers 
placed at the common mid-point. 

The purpose , is to, add records .illuminating, the same 
5 - .point in the. subsoil s to increase - the^signal to,:--npise 
ratio, and,, the , primary reflections : to . secondary 
reflections, ratio, thus benefit from'; the advantages ' 

of "multiple coverage". '•>•'■"■■, 

• In ° rder to ™ake ; . a ' zerb ^offset image, a so-called 
3 oblxg^eness ; . correction method^ pr £S (Normal Move, Out) 

.. : correction method is used. , 

If it. is. assumed that the subsoil is horizontally 
stratified with no, lateral variation. of f .the propagation 
velocities, it can be. shown that records that have the 
property of. illuminating the same point in the subsoil 
are records that, have the same .midpoint. . , 

. How ^. r .' the fmage of a reflection in. the subsoil 
. arrives at , variable^ times depending . on the offset 
Therefore,. ln order to stack reflections, it is necessary • 
to start by correct^ tne .arious; record, to Bri ng \t h e m 
all to a common 'zero offset reference., , 

, Historically, the NMO correction is phased on a 
■ particularly simple model 'of " 'the "subsoil ; ' ^homogeneous 
model with horizontal reflectors. ..... "' f ' 7:-T v " 

f , In AfUs,. J^del ^ ; reflections^ associated with " each ' 
, subsoil reflector ..are theoretically.; aligned . along 
hyperbolas also called indicators, centred 'vertically 
above the midpoint/ . A ^y\ 

The arrival time of a reflection is then a 
hyperbolic function of the source-receiver offset, the 
•shortest time being obtained at zero offset. 
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15 



in whxch ^represents the^zero offset travel Wme; and V 
denotes the average propagation velocity of waves in the 
subsoil ."' ' *.•«-■• 

The simplistic model mentioned above is ? based 
particularly on assumptions such as low angles ' of 
incidence and an isotropic medium. 

But assumptions ^made fdr' the simplistic ; model are 
, tod restrictive to ' describe a complex medium, and in 
particular cannot be applied to the propagation of 
^„*yb*W in v in anisotropic 1 Wium n Mmedium.:ih; ; which - - 

the : 



2 0 propagation 



The u S e' of : 'a * lVs V sinfciistid model " ' is- necessary, '. 

• particularly due'to: ' ...... ; : 

' "£ the:, use of. -long ^streamers for acquisition iii; deep ■ 
offshore., which in : particular means that traces with far : 

25 of f set s ' are acquired ; / 

L observation ' of 'anisotropy in clayey type 

sediments . 

The generally . accepted assumption consists of 

modelling an anisotropic medium as a stack of 

30 transversely isotropic layers with a vertical axis of 
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in order to . stack records in each gather, the NMO 
correction straightens the hyperbolas until they are ^ 
theoretically horizontal. 

The- NMO correction is then made based 'on the 
following hyperbolic ^equation- for the travel ; time Rafter 
reflection,^ assbcia^ Source-receiver of f-set- x 

. pair':" u/ . * ' " V4, "N" % *- ^ / 

X 
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symmetry. This is known as Vertical Transverse Isotropy 
(VTI) . • . .. 

Thus, it has been proposed to determine NMO 
corrections to be made by: • ^ • 

5 " ' : " introducing vertical non- homogeneity r;i'ntp ■■ <"■ a 
homogenous VTI medium model, V as presented in the document 
"ALKHALIFAH , T . and TSVANKIN I., 1995, Velocity analysis 
for transversely isotropic media... Geophysics, 60, 1550- 
' 15 66"; ' . '. -"' ■• ■- ■ ''. - v ]\ 

"* or by " ' * v. ^ ' , « 

.-, - introducing VTI anisotropy . into a model with 
stratified isotropic layers, as demonstrated in document 
" SILIQI R - , and BOUSQUIE N. , 2000, Anelliptic time 
processing based on an offset hyperbola approach, 70 th 
Ann. Internat. Mtg. : Soc. Of Expl . Geophys ., 2245-2248 » . 

This second approach combining vertical non- 
homogeneity and VTI anisotropy to give a new model of the 

^° llf aPP ^ S 'M-^W^h^ ? ^PP rb ach ; .in.;;'most: real ? 
- cases studied.- = ; ; .', ' '. '• ,! " • 

>" ] The following equation', for correction "of the "travel 
time, t; after reflection/ %sing"' an ' anelliptic offset ' 
hyperbola, 'associated with 'a source' ^ receiver ' of f set 
'• . pair,.' is derived from this model: • ' ' " 

where, V. is the velocity^ convent ibnally used in seismic 
corresponding to small offsets, and n is a. parameter, 
called the anellipticity parameter. 

The document "SUAUDEAU E. and SILIQI R. , ,2001, 
Anelliptic pre-stack time migration, Annual International 
Meeting, CSEG Expanded Abstracts" also proposes to 
include the NMO correction by anelliptic offset hyperbola 
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in the equation for the path time used during the PSTM 
migration operation. 

The PSTM migration equation - is conventionally 
expressed in the form of a double anelliptic offset 
5 hyperbolay . .which is the sum of two square roots (the 
Double Square Root (DSQR) equation) .. - iVV^.-'> 
, : . Taking, account:" of an4iiipticity, ; } the expression for 
this constant offset equation, becomes : *, - ''■?*■ ? - f : 

8T1 t , If tn/Q 2 - (X-X.Th) 7 (fWlY^^Iz^ 

t = 1 + 8^ vtl+irjj (l+8ri)V* tl + Zn): 0 + 8n)V* 

10 - : " \i : '•: "Equation (lb) *>•>..■ 

where : 

- the V and r\ parameters are those mentioned above, 

- Xn is the coordinates of the midpoints, 

■ - x - x™ represents the migration aperture, 
15 - h . is the half source -receiver offset, 

- t 0 is the double, time, for zero aperture of the 

• •. ,' operator. . ' . . r -y _.?■..;>...•'.. y-_i J.- ' 5. --.v-"-. >. •- 

V. Note that when ; the aperture x - x™ of , the migration 
/• ' iV^zeroV- the ''PSTM' correction equation . (lb) becomes., .the 
20 NMo' correction " equation * (la) . Thus, the NMO correction 
forms a special case of t the PSTM migration; the case of 
PSTM migration with zero aperture. 

Therefore, . .finally,! in order, to: perform seismic data 
processing including^ NMO correction taking account of 
25 vertical heterogeneity and VTf type anisotropy, it is 
necessary to determine the velocity ;v and aneilipticity r\ 
parameters . 

The said parameters V and n can conventionally be 
estimated by carrying out two passes such that: 
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- during the first pass, the distribution,:, of 
velocities V along the . time axis is estimated/ only using 
data with near offsets; 

- during the: second . pass f - the anellipticity r\ is 
5 v estimated along [ they:t ime : Vaxi : s-y 4 lis ing •: ^.v .,>;-v^i... • 

the ■ distribution of y velocities ^ : ^deterniined 
during the first pass; ank^fr v . Vr> v.. 

- all data (including^ data^ with far^ offsets) " 
. The, document- entitled: oi^V JfLAv ^-i * 
10 - "SILIQIa R:L ; 20 ;; P1, $s Technological^ : leapt* f in-,,; time 

processing focuses the data throughout- anisotropic media:. 
First Break) .19, No. .11, 612-6185 a U> « - \ ^ 

has also demonstrated that the parameters V and T] can be 
estimated in a single pass,' during which' bispectral 
15 analyses are carried out in order to simultaneously pick 
the V and tj parameters ; along the^ ; t ime axis, using all 
• •.. • . v data\-\- ; V->*'^ :t:^^v.;.^^-.^c4 . *'~& : \ 

However, dense analyses of NMO' correction parameters 
are 1 ' preferably-; made' when- the - correction equation • no 
20 longer depends : Toh'" the,. time^t 0 ' -(we' v their talk about static 
NMO correction)'- * ' ' - - - - , - - - > 

A static correction effectively pr6vides a^ means" of 
- offsetting all sample s^eaclf forming traces £' -of tHe same 
.: time 8t, for a given' offset 1 ; u ;.V« \ •fr7^^r:y-^; / ; l; • \ , 
25 Thus, when : a static correctidrt is made, the- number 

of calculations- to be made"; can 5 be signif icaritly reduced 
and the trace Stretching* "phenbtheuo^ eliminated, which 
makes the said dense analyses viable. <''-.' 

At the moment, the only way to make velocity and 
30 anellipticity analyses is to use two passes, in 
particular making use of parabolic approxi mat ions of NMO 
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residuals, in order to obtain a dense picking of the V 
and r| parameters. 

In this context : 

- the : first step is to estimate residual velocities, 
using near offset data ,v according 1. to= a first estimate ; of 



velocities; w • i : W ••. ? " *i 



■v 



- the next step . is . to estimate . anellipticity; on all 
dat a , • usin'gr velocity^ updates* done previously, t \ . 

A mute .law must also be def ined: ; to;^stimate "residual 
10 velocities; so asi to ' v only keep data that are . considered 
* as being near of fset data . ,: ; *<; f -; v « i -r-.v vv " w . ,. §v 

However, the efficiency of the two-pass'' analysis is 
particularly, sensitive to the choice of such . a mute, law. 

On the other hand, the anellipticity r\ is estimated 
15 mainly making use of far offset data. 

, c: But, the « anellipticity --est imate >...made.' for .far offset 
data is not very precise, such 'that the correction made 



is finally^ inaccurate;-^ V 



* ^ One- purpose; of t ^the^ invention is to overcome these 
20 limitations and disadvantages ,,- by : proposing^ a : . processing 
•.•including a static NMO correction which is^ more/ : efficient 
'and more * precise \ ^ ^ - . . ' T j : ' . 

... .More, precisely,^ the purpose of; the invention is. a 
dense ; determination: p ' aneil ipt i c i ty T| ' 

25 parameters in a single pass % usingf all available data, in 
other words, making^ use of all, offset ranges r . , 

, Tq ; , achieve.?' this , : \ the invention proposes a method of 
determining the velocity V and anellipticity r\ parameters 
for processing seismic traces from a - common midpoint 
3 0 (CMP) gather including an anelliptic NMO correction, 
comprising: ■ 
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■ - a preliminary step to define a plurality of nodes 
(dtn, t 0 ) , the' said nodes being indicative of parameters 
dtn and x 0 representing the NMO correction for the 
maximum offset and the zero offset travel time in 
. 5 hyperbolic coordinates,, the . said, preliminary step being 
followed bv . . . : v ... 

: , for .each node,, (dtn, r^o)^ defined, in the preliminary 
step, , steps for the NMO correction of traces ' in the : CMP 
gather as a function , of. , the .values.;. : of the said, parameters 
10 dtn . and, x 0 . at the,, .node considered; and f orv calculating 
the : semblance function associated:- with , the said NMO 
correction for the node, considered; , and ^ 

- for each picked time t 0/ a step including 
determination of the maximum semblance node (dtn(t 0 ), 

15 T 0 (t 0 ) ) , 

- and a final step to convert the dtn(t 0 ) and x 0 (t 0 ) 
parameters so as to obtain the . velocity^ V (t 0 ) ; -:. and 

>... anellepticity .r| (t 0 ) laws.'^v,'^ , : .1 . ?v: v r ;,. v ' . 

: \^ v According ^ to a .first* embodiment* of t the invention, 
20; the processing ;done, ; :is . a r - • & t at ic f .NMO. ^correction • -ofe . se i smi c 
" ■ traces 

According to a second embodiment of the invention, 
the processing done is a PSTM migration of seismic 
traces, the said, PSTM migration including a static PSTM 
25 NMO correction of the said seismic traces.; 

A preferred but non- limitative aspect of the process 
according to the invention relates to the definition of 
the dtn and t 0 parameters relative to the velocity V and 
the anellipticity r\ so as to make static NMO corrections, 
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according to dtn = l — t Q +J 0 + — and t 0 = 



l + 8r| ' * • J ; 7 : 

Other aspects ^ r .purposes arid advantages of this 
invention will become clearer after reading- the following 
5 detailed description with / reference, tortile "attached 
' Figures in which: • - y '"* U' vV " "* > ^""^ — • 

- - Figure la shows the anelliptic of f set :: hyperbola 
used to make the NMO s ; correction : and" illustrates^, the 
meaning of the parameters t 0 and dtn; r \ : 
10 - Figure lb represents the DSQR equation of the 

double anelliptic offset hyperbola of the PSTM migration 
and illustrates the meaning of the x 0 and dtn parameters; 

- Figure 2 illustrates the effect of the t 0 
parameter of the reflection curve corrected by the NMO 

15 correction/ ,. , v. 

- Figures 3a a.nd 3b represent the analysis volume 
; (t 0 /'dtn, T 0 ) in which the bi spectral pi eking of " the dtn 

.- and x 0 ; ' parameters' according to the invention^ is\ made; .V 
• Figure 4 combines each of (V, Van) and (dtn; x 0 ) 
20 approaches, representing their corresponding bispectral 
analysis panel; 

- Figure 5, shows the correspondence between the 
(dtn; x 0 ) and (V, Van) ; pairs of NM6 correction parameters ; 

- Figure 6 shows a CMP gather of true seismic traces 
2 5 before the NMO correction and the bispectral picking of 

the (dtn, t 0 ) parameters corresponding to this gather of 
traces; 
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- Figure 7 shows the semblance function and velocity 
and anellipticity functions deduced from the picking of 
the dtn and x 0 parameters in Figure 6 . . 

- Figure 8a is a' flowchart representing the steps of 
5 ' a " first' 'particular ^ertopdiment off' the' invent iox?7"' namely 

'determination 'of the "v and" rf parameters to make an NMO 
correction; >: " : ' T*'* ' 

^ ^ ^ Figure sS is a flowchart repreilntlng^the'steps of 
a second particular embodiment of the ;: ' r invent ion namely 
10 determination* of^ V: and r| ' parameters to do - a PSTM 
migration; : - V :> ■ * ' v •• y ' - 

- Figure 9 is a diagram illustrating the different 
~ 2 operations carried out in order to determine the V and r\ 

parameters to make a PSTM migration. 

15 In general, the method according to the invention is 

a method for processing r records ? of variable offset 
seismic traces, and- this processing uses the recorded 
seismic traces to build common midpoint (CMP) trace 
gathers, arid the - tr^ce^ in each; gather^ - are subjected to 

20^ an NMO correction.^ . ^ ^ > * ^ r - 

; In particular , the . method according, to„ the invention" 

. determines., the^ velocity; V_and; anellipticity: rj "parameters. 
. to perform processing including 
the seismic traces of a CMP, gather. ,.. \ 

25 For example, the said- processing may be: 

- an aneliiptic NMO correction of the seismic 
traces'; • F — ■." • 

- an aneliiptic PSTM migration which, as we have 
seen above, jointly applies TZO and migration operations 

30 (in the following we will refer to a PSTM NMO 
correction) . 
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The following- description is more specifically 
applicable to the NMO correction. However, considering in 
particular the two particular . . . embodiments of the 
invention that will . be /.described, below, we can see that 
5 this description is. equally;: applicable to any. processing 
including an NMO t; correction, particularly processing 
including a PSTM NMO correction. * . * "; 

/Two new . parameters .. dtn, and - i 0 : are, considered in 
order to, specif ically determine V : . and r\ : . • .• -. w .. 

10 : : ' cv'c ~ T o. .that; represents; the zero offset travel .time h in 
"hyperbolic coordinates" ,( see Figure 1), '/^Ky. 

t 0 = — — r[ Equation (2) 

1+8-n ..; : 

and, 

- dtn that represents the NMO correction for the 
15 largest off set. (see Figure 1) , 

dtn = t x =xo« ....t t x= p namely: 



dtn ,= -^-it 0 +'. 



U + 8ti^ (1 + 



max^ Equation;. (3) 



• ,- : It;,is important • torncrte- that: dtn: isVdef inei relative 
to the- - velocity V and the ane 1 1 ipt i c i ty Y| , while x 0 is .. a 
20 perfectly ahellxptical : parameter defined relative to the 



anellipticity r\f independent, of V. SK 



Accordingly, the 4 velocity V : and - anellipticity r\ 
parameters may be calculated. in accordance with equations 
(2 )• arid ; (3) - using 'the: following 'conyiersiori Equations (4) 
25 and (5) : 



V = Xrmx , Equation (4) 

\dtn(dtn + 2r 0 )^ 
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and ii = -M — -ll Equation (5) 



Figure la shows, the . jan^^ hyperbola 
. 5 • used to make the - : : mo" correction and . illustrates ' the 
meaning of the x 0 : and -dtn parameters.. : . 

The said "hyperbolic^ coordinates H are shown xn this 
Figure 1. Their origin on the time axis is taken, at the 
intersection of . the said time, . axis with . the. . asymptote . ; 
10 tangent to the said far offset offset hyperbola (see 
equation (la)) . 

Using the (dtn, % 0 ) parameters,, equation (la) for 
the anelliptical offset hyperbola becomes: 



J 2 dtn(dtn + 2x 0 ) 2 „ ^ , , _ x 
Tq -f — — t^x Equation (6a) 
• X max . 

15 The (dtn, t 0 ) parameters defined for the velocity V 

and ;M anellipticity r\ are thus used ' to . make , the NMO 
•? i; correction CORRnmq: = t- — : 't 0 " to be" "applied to"' of f set yx< 



trades "independent of t 0 : 



GORR^U) ^ -to' V xg + 4tD(dta+ , 2to) x 2 Equation (7a)^ - " - - 



20 .. .-• Therefore, this, is a static /Vl^O^cbrrecti.on'^ritt other 
words, the data recorded ori a given offset trace will all 

.. ' •>.< ■• «••.-. . . .... . . {?/..* ,v.s : ; Viv*'- : •• •' ••'■^ , 3'" ""' : - c ':-'-. '::*> v : •'".•v , . v -'' ":f ' •/ • ' > ; -. •- 

.' be corrected in the, same . way A f orv (dtn, , x 0 ) pair, 
independently of the time at which, these data were 
acquired. 

25 Consequently, the estimate of the velocity and 

anellipticity parameters is not disturbed by stretching 
of traces usually observed when dynamic NMO corrections 
are made . 
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Similarly, using " the two parameters (dtn, x 0 ) , 
equation (lb) for the DSQR double • of f set hyperbola for 
the anelliptic PSTM migration becomes: 



;jlo : , s dta(dtn fgtoXx ^x m Vh)* 



t = t 0 - X 0 + 



• 1 - : >: :v ?> ■ feti ^fP ~ h Eguatapn (6b) 

•*. • max--.:/-. ■ 

: •• ; - . • • i • - i& . .:.f v * • .> t -o'-,:4^ ^ • . . • ' 

Figure lb shows the DSQR. double offset hyperbola of 

the PSTM migration defined with the (dtn, x 0 ) parameters.. 
10 In the context of the PSTM migration, x max represents 

the maximum offset and aperture of the migration. . 

The (dtn, x 0 ) parameters defined relative to the 

velocity V and anellipticity r| are thus used to make the 

PSTM NMO correction CORRpsm = t - . x 0 to be applied to the . 
15 offset x traces independent of x 0 (equation 7b.)>: ; : 



20 



.Gorr 



s?m(x% = ^ EE 4to{dtn + 2t 0 )(* -x- hf\ . V 

V 4 : ?T ] V 4 • • : * x 2 * \ " -~ 

Therefore, this is a static PSTM correction/ In 
other words,, the se of aperture samples x/'-' ;;x^ : »of 1 an 
" iso-of f set" cube are offset by the same time for a given 
(dtn, x 0 )' pair/"'" ^'r V/ ; ^' v< ~ 7 r f • :> ' 4 * 

By analysing the f (dtn,;- x 0 ) V parameters m several 
picking times/ the method^ to the invention can 

be used -in particular to T determine the ' (V, rj) parameters 
necessary for processing including an anelliptic NMO 
25 correction of the traces of a CMP gather. 

The . said process includes the steps presented below 
in a simplified manner. 
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. During , the preliminary step, vi an analysis volume is 
defined including several nodes (dtn, t 0 ) - 

The following are performed, for all nodes in this 



volume : 




Finally, the following are determined for each time 
in the several picking times: 



15 - firstly, the node (dtn; t 0 ) used to make an 

optimum correction, for example with regard to the 
semblance criterion (this type of criterion typically 
being used 'in ' seismic/""" processing' to "measure 1 the 
horizontality" * of reflection* curves- and to determine the . 
20 reliability of the : picking) / ~ -'i > >; - , -V v ^ 

- secondly, the values of the dtn, t 0 parameters at ' 
the saidi maximum semblance i:node/ aret converted to'" values 
of the velocity V and anellipticity, r| parameters at the 
'. said picking, time' ponsidered., ^ : v ; t / * v v ..... , . 

25 The velocity law ; ( in other words all (picking time f , 

.V)) pairs and the ^elliptiicity f law ... (all,, (picking time , 
r\)) pairs are thus; defined ;v .i : 

Finally, the NMO correction for all seismic traces 
may be made using, these velocity V and anellipticity tj 
30 laws in equation (la) for the anelliptic offset 
hyperbola. 
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Similarly, and as will be described in detail later, 
the (dtn/ t 0 ) parameters of the static PSTM migration 
(see equation 7b) ) may also be checked: The velocity V 
and anellipticity n ..laws are then determined and may be 
5 used to make ''the i^TO'^migra'tion^ 'l^the^DSQR equation (lb) 
: ' f or the double anellipt i c ' of f set ''hype'rb'ol a . . ; v > 
• " Figure 2 illustrates the effect 'of the parameter To ,. 
(and therefore the" anellipticity ^ "according to equation 
(2) on. the curvature residues after the NMO correction. 
10 ' Note that the vertical . time scale, of the curve in. 

Figure 2 is' exaggerated so that this effect can be. 
clearly understood. ■ v 

Three curves are shown in Figure 2, for which the 
dtn parameter is fixed to the correct value and the 
15 parameter To is assigned different values. 

The central curve shows . the case in l which To is 
/-equal J to its f correct w value Toi .in ^ other, words ^ when ,.0ev ' 
. ^ ^corresponding;, anellipticity T)x .is, equal ^o/ its ; true, value.- v 
Ti tru e. As expected, the' corrected; reflection" curve.. is, then < 
20 ; horizontal.' • _ :* - . . >" . ' - ( .V . " 

' " - • The upper curve " shows /the . case in which t 0 ; is? equal 
to a value t 02 less, than its . correct, value t 0 i, the 
corresponding-^ 
: true value. 

25 The ' lower" "curve represents the :: case : in which t 0 ' is , 

equal to a value t 03 greater than its true value t 01 , the 
anellipticity n 3 being less than its true value Titrue- 

Note that these lower and upper curves show that the 
"horizontality" of the corrected reflection curve is 
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acceptable at. near offset (x « 0) and at far offset (x « 

On the other hand, significant curvature residues 
are observed when, the offset x no longer tends towards 
5 one of these limiting values. t 0 and^ax'. In particular, 
particularly significant residues are observed for an 
offset x centred in the middle of the offsets range. 

For example, and as' is shown in" Figure 2, when t 0 is 
equal to x Q2t an RMOt 02 residual correction must be made 
10 to the offset traces XrmoT 0 2 • Similarly/ when x 0 is equal 
to x 03/ a residual ' correction Rm6t 03 must be made to 



offset traces x^c/toa • 

Due to these significant curvature residues, the 

offsets range can be used almost in its entirety in order 
15 to determine the aneilipticity r[. 

Therefore/ setting parameters in (dtn, t 0 ) for the 
r NMO correction make it possible to use available data for 

all offsets (x between 0. and x^ax) . in determining^ the 

aneilipticity tj. 
2 0 •* As. has already been xnentionedv • this : is not the case 

for NMCT corrections' jEorc which parameters have- been set 

with V and r\ for which the estimate of the aneilipticity 

r| is made essentially using- far offset -data . v. 

. Thus,- the effect- of t he V ne wV ane 1 1 ip t i c pair ame t e r To ' 
25 is distributed on all of f set^, ; imlike the aneilipticity i\ 

that only affects fair offsets . Therefore , the "behaviour" 

of the parameter t 0 gives a better constraint on 

aneilipticity values. 
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As already mentioned above, the optimum dtn and x 0 
parameters are determined within a 3D analysis volume 

(t 0 , dtn, T 0 ) . 

Several nodes (dtn, x Q ) , . in other . words several 
5 pairs of dtn, ..t 0 ' parameters for : which ..values, are known, 
are considered in the, said analysis volume . ^ :j . 

Nodes are usually ' regularly., spaced . from . each . other,,, 
by an increment ' Adth on the dtn axis and an increment At 0 
on the t 0 /axis v • '.-.J: - • : : " V- /-; : . - ,/; .-/- ; \- 

10 Minimum values dtn m in, T 0 min/ t 0m in and . maximum values 

dtn max/ Tomax, t 0 max of the* t 0 and t 0 respectively provide a 
means of defining the limits of the said analysis volume. 

Advantageously, plausible values of the velocity 
parameter V and the anellipticity parameter i] may be used 



15 to define a corridor .[dtn m in(to) , dtn max (t 0 ) ] , [T 0 min(t 0 ) / 




Figures 3a and 3b represent the said analysis volume 



(t 0 / dtn, T 0 ) in the context of an example of true seismic 
data processing during an NMO correction done according 
3 0 to the process according to the invention. 
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Figure. 3a represents three 2D panels a, b and c of 
the analysis volume: 

- panel a is a panel., (dtn, x 0 ) with constant t 0 ; 

- panel b is a panel (dtn, t 0 ) with constant x 0 ; - 
5 - panel c is a panel (x 0 , t 0 ) with constant dtn; 

..Figure 3b diagrammatically s : . shbws'* the 3 analysis 
volume (t 0 , dtn, \ x 0 ) ' aiid.. three * intersections with this 
A volume along three planes with constant t 0 , : x 0/ dtn 
respectively, eacii \of : -'these' intersections being projected 
10 ' onto the corresponding side a /• b or ' c iri Figure^ 3a . : i f : . ■ 
Panels (dtn, x 0 ) with constant t 0 (panel 1 a in the 
above example) are 1 panbls' in which the bispectraT picking 
of the dtn and x 0 "parameters is made, for Example 
according to the maximum NMO correction semblance 
15 criterion, for the considered picking time t 0 . 

The [dtXW \t 0 T, n dtllniax (t 0 ) ] " r [X 0 min (t 0 ) Xomak (to) 3 



corridor mentioned above" for the" effective" analysis 
■ ' ■y inside , the analysis^ volume ds ^ also - shown " in Figure d 3a 
(lighter sorieon each panel)". v\ r * ~ , i ,v „.* 4 

2 0 : , The s e i smi c ; dat a samp ling t ime s t ep A to r de f ine s the 
difference between two ' successive 1 picking 5 ■ times ? (times 
for which the maximum associated semblance node (dtn,' 5 t 0 ) 
is determined) and therefore the number of bispectral 
picking, panels to lie considered.' ^ K. 



■ - * 



.25 .: Automatic ' picking^ also provides a means of 

. extracting the " dtn, f 0 parameters at a ' density that is 
greater when the At 0 increment between picking times is 
small. 

Sampling of the analysis parameters dtn and x 0 is 
3 0 directly related to the resolution of the . seismic 
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exploration; dtn and x 0 effectively have the same 
dimensions as seismic records . 

The • systematic search for the maximum semblance, 
conventionally known in itself, provides ..a means of 
5 determining the • ; <dtn,>; x 0 ) ; pair providing best focus, for 
a given picking time to . • - , v. ■ ^ ^. - 

Parabolic inte^olations about the ...values, of nodes 
(dtn, / x 0 ) can also provide a means of evaluating values , 
/.of the dtn, t 0 parameters • between-:^ the., dif ferent , nodes 
10 that have - actually v beeni;pipk^dU; Andv in^vparticular ...this 
type of evaluation makes determination even more precise 
. (in contrast., to determination limited, to nodes, in the 
corridor) of the dtn, x 0 parameter pair maximising the 
semblance function. 
15 Finally, the velocity .V and.,- anellipticity r\ 

parameters are determined, always for the picking time t 0 
considered, .by using.' the* above, mentioned conversion 
. equations (4) and : ^(5). ; - « ~ * >\\ : 

Figure 4 shows : the ; comparison,; • between . t^e^ , two 

20 approaches" (V, Van). :< and ,(dtn, ^T 0 ) , by representing the 
bispectral picking panel, for each, for a given picking 



•time.. v <=■.-. . -., 



The figure at the right illustrates the ^.conventional 
approach (V, Van) for which ^t he . two axes . are the. velocity 
25 axes (the anellipticity r\ being, related to the ratio of 



1(V* " 



these two velocities according to T) = ^l"^" 1 



) • 



The left figure illustrates the approach (dtn, t 0 ) 
according to the invention for which the two axes are 
time axes- 
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It is. important . to note from the study in Figure 4 
that the dtn and x 0 parameters appear to. be decorrelated. 
This "decorrelation" . is. striking when the two approaches 
; are., compared,- spreading . of the :. ,(dt^. Jp)r-: ; spectrum, 
5 actually being much ;.^rr^ 

v Van) spectrum.- : ' • * f ' 

Therefore the picking made., in. the context pf ' ..the 

• (dtn, t 0 ) approach, according , to^ the \ invention- is more 
precise than the „ approach., .. 3 ca^r led .^put in gon^ntional 
10 methods. 

\.. Furthermore, this decorrelation enables filtering of 
the dtn and Ta pickings separately, while keeping NMO 
corrections. This is, not the? case for the V and. T| 
parameters for which a reduction in one of the parameters 

15 must necessarily be compensated by an increase in the 
other and vice versa.. 5 , .< ■* . * * 

:'.. } . And due to interpolation^ 
of ^parameters according to: ; the^;inyention /V dtn-yand^^- it. 
is'/the'ii; possible to* perform- simultaneous - interpolation 

2 0 and f iltering of the standard, NMO. • correction parameters V r 

and T|. ; . 4 - - V 

Figure 5 shows non-linear correspondence:^ according 
to equations (4) andy (5) mentioned" above,- between the 
pair of time parameters • (dtn, T 0 )i' ; and , the ; pair \ : 'of 
25 velocity parameters , {V> 4 Van) r * , 

. Figure .\ 6 . shows.;; fforn left- to} right; :.; \ ;V ; l \vi : :^>V"v 

- a CMP gather, of real, seismic * traces^ before -NMO 
correction; 

- the picking of the dtn parameter corresponding to 
30 this gather of traces; 
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- the picking . of the x 0 parameter corresponding to 
this gather of traces. 

The straight line x 0 = t 0 * on the picking of the 
parameter t 0 at; the right in Figure 6;*- corresponds to the 
5 purely hyperbolic reflection- curves • J ; ' : 

From right to left/ Figure 7 represents the 
semblance function \ and*" the 4 velocity - V and the 
* "anellipticity r\ functions • (see equation (4) and (5) ) 
deduced from the pi eking* " of ' the dtn and x 0 parameters 
10 shown in Figure 6. * ..^-m 

On this example of real seismic data processing, it 
can be seen that the values 'of V and " T| obtained generally 
corresponding to a semblance of more than 40%. 

The following description contains details of two 
15 particular embodiments of the invention. 

The ..first of these modes relates v to a process for 
determining optimum parameters to make an ahelliptic NMO 
/ ^ correct ion// to " the traces^ of , a</ CMP - gathers: (see the 

' diff erent " steps'; shown - in-- the; flowchart-' irr Figure 8aK 
20* w With reference to Figure 8aV this ' first embodiment 
includes an initialisation- step* la - during which r the 
following; operations^ are /carried ■/■•out in; sequence:; v\ 
> ,:' -\ determina ; tioiil. of the - limits of the analysis, volume 

•;/"v- [ d t n m in-i . : d t n m ax ] :V. • . [ T Omin .•/ • t Omax 3 V; . [ t omin / 5 t omax 1 J S'tv- >:'^.v ; s. ; -,:i ■ 

25 r calculation of NMO corrections CORRnmo (equation 

... (7a)) for all offsets' and, for all,, nodes (dtn, x D ) 
included iri the analysis volume; ,,T 

- delimitation inside the analysis volume of the 

corridor [dtn min (t 0 ) / dtn^ tto)] [Tomin (t 0 ) x 0ra ax (t 0 ) ] of 
3 0 plausible velocity and anellipticity values. 
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Once the initialisation 'step la has been carried 
out, a step 2a to calculate velocity V(t 0 ) and 
anellipticity T)(t 0 ) laws is performed for each " gather • of 
CMP traces. : ""' ;T ' 1 " ' • .-^ " " '* ' /' !' 

5 1; ' - This 'step' 2a ^cornpfises :^ I ^ ^:/^;.^f ? ;v : >• v 

- a first operation^ 3a;| carried 1 " out At f or each node 
(dtn, x 0 ) in the corridor'def ined 1 in^ the ^ initialisation 
step la; ..during :v -which ?> the v following -** operations are 
carried out' in sequence for'; each picking" time to : * " " v 
10 - application of . static NMO corrections CORRiJmo 

for* all offsets along the corridor? precal ciliated during 
the initialisation" step la; - ;> 'i^r<^% : * ii. 

- calculation of the nsembiance^ f unction on ; data 
corrected along the corridor using a time window 

15 appropriate to the dominant wavelet; ; 

- summation ' ( "stack f,v calculation) • of data 
corrected along the corridor (only near offset .data can 
be used advantageously for this^ purpose) ; 1 %r - 

• a second: operatic 

20. time, t 0 „ (the said, times \ being, at intervals of . At 0 between 
- ■ [tomin/j tomax] during,,, which . the following operations are 
carried out: " 1 * • - ' <: / - • . V'*; 

- search for, maximum semblance in v the; [dtrimin (t 0 ) / 
, dtnmaxito) ] [xomin (t 0 ) . Xbmaiai (t 0 ) I t corridor^ ofer the, bispectral 

25 paiiel- (dtri, t 0 ) ; - v ^ * ^ ' • 

- check the fact that; the position in (dtn, t 0 ) at 
• the maximum semblance corresponds to a summation extreme 

value for the same values dtn and t 0 ; 

- creation of the dtn(t 0 ), to(t 0 ) and semblance 
30 ' (t 0 ) series; 
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- a thii-d operation 5a designed to select and adjust 
the pickings- obtained, during which the following 
operations are carried put : :> _ 

- increasing sort of the semblance series ..(to).;,... 
5 validation 'bf pickings^* dtn. : ^ 

.... distance in time' tq ># the highest _ j . semblance pickings is 
* greater, than^ a, predefined, value; ^ ? .,. : ^ , 

'TU, " a&justn^^ x 0 : values 

validated by parabolic interpolations A . using surrounding 

10. , values; ... ^w-,^ -^m Xs?\:t.,* -• ,- 4^ * - . 

• - retention, of * picked, , validated and adjusted 

values if it is possible to , calculate ^ the. Dix interval 
velocities- -.;■■?.; with .-■••the pickings .*' ; with -the highest 
similarities. . ■■■■■ ' . 

15 - a fourth operation ,6a designed to use : equations 

: (3) and ( 4 ) to cpnyer t, c pi eked , validated and A adjusted and 
c-: retained ^values : o.f x:; dti& and Tqvdur ing-t- operation:;-.,;-. 5a>.r :into . 
velocity V and anellipticity r\ • laws . ' ;; > • ■ ' "^ A * 
' . : " ; \ V . The" velocity V- ' and" anellipticity^ as / . a, 

20:: ; furictiori'^ of time -are thus- perfectly^ - determined The 
• — « anell Iptical^ NMtf'^ correction— of "seismic traces < in the CMP. 
' ' gather;, caii'' thiis r Be done-'- preci sely . f : ;*\ ' " ' v *' 

Th^^ second^ particular embodiment:; of v the invention 
relates^ tox. a-; process ■. for determination of optimum 
25 parameters for anelliptical PSTM migration pf traces^ in a 5 
. ' CMP gather. . ( ' A \, . 4 \. c , / 4 \ ' 

This second , embodiment may be included in a 
generalisation of the first mode discussed above. 

As has already been demonstrated, the use of the 
30 (dtn, t 0 ) parameters, enables static PSTM NMO corrections 
(see equation (7b)). 
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In the context of a PSTM * NMO correction, 'this 
application has the same advantages as previously 
discussed for the NMO obliqueness correction. 

More precisely, it will be noted that the first 
embodiment is.^ only a special ' v case^bf * * the second 
embodiment i corresponding to the case 'of a "zero migration 
aperture. . ' • .< v • • ;H ;r „^' > Vt ; V 

With referencing 
includes an initialisation .step lb " during which the 
lowing, operations are carried out in sequence- tv ^ ; ■ 

- determine, the. | limits- of% the f analysis volume 

[dtn min/ dtn max ] , [Xomia, Torna*]/ [t 0 min, t 0T 4x] / " * 

- calculate NMO corrections^ CORRp STM ^ (Equation (7b) ) 
for all nodes (dtn; x 0 ) ' included in the analysis volume 

15 and for all migration offsets inside the migration 

- delimitation of plausible v " velocity " and 
ahe.llipticity values^ inside' the ^analysis Vofiime^ of I the 

• .corridor [dtn^(t 0 ), dtn max (t 0 )V fWn (t 0 ) W " (tf 0 ) M ^ 



10 '.: f0ll< 



?°. ' Once t His initialisation step- lb ' has been" done /' "the 

■ ' said firsts' .embodinjent -creates a' : step _ to' calculate the 
velocity ' V(toh;and anellipticity r\ (t 0 ) v laws' >or^ 'each 
gather- of: CMP traces x~ ' . * 1 - : ° ov: AoX'* ^ V 

.This step 2b comprises.: ,r - *'~ " . >< 
25 - a first operation 3b 'done'' foif each node'*'' (dtn, x 0 ) 

of the corridor defined ; during the ;; ihitialisatioh step 
la/., during : ". whi ch ; t he foil owing are 1 carried * but in 
sequence : '"""■ ' 

- for each offset class, operations for : 
30 - application on all midpoints inside the 

migration aperture along the corridor, of static 
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corrections CORR PS tm - precalculated during the 
initialisation step lb; 

- summation of corrected midpoints along the 
corridor. ... ,, jf J./.-;* ■; . 

5 ■., , . r for each .picking operations to: ' ;r 

v; ...^;>y'. lV< . n.calculcl^ f unction " : on 

corrected data along the ^ corridor using a tim^ window 
appropriate to the f dominant; wayelet; •.' ' • ; 

r - summation (stack . calculation)*. of < data 

10 corrected, : along . , the,^ corridor v>(only data v with small 
offsets, can be. used /advantageously .for this purpose) : 

- a second operation 4b carried out for each picking 
time to,, similar, to , the^ operaticm 4a described; above, to 
create the dtn (t p ) , x 0 (t 0 ), and semblance t 0 series; 

15 - a third ... operation similar w to operation ■ ■ 5b 

described above, designed to select and adjust the 
pickings obtained; • 

- a., fourth ..operation/ 6b similar to operations 5a 
described-, above;, designed to, convert .Lvalues^ ofo dtn/ x 0 

20 into velocity, V, and anellipticity tj laws. . . ^ 
. •' <! . The,;,- velocity ' v/ and ariellipticrty -ti 1 raws ^ as ^"a 
function of- time are thus*, perfectly, determined /* And the 
PSTM migration of seismic, traces - in the CMP ga their • may 
• ' thus be ^ done accurately . y^\ 1 / « ^--'^ « ; rjf 
25 Figure 9 i 1 ius t f a t es, ^ the second embodiment* of the 

invention that has just been. described r • ' T 

Seismic data are initially grouped /into *. isp-offset 
cubes. 

As has just been described, NMO corrections CORR PSTM 
30 are applied (see first operation in step 3b) for each 
offset class (namely for each iso-offset cube) on all 
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midpoints . The arrow marked with the MIG label in Figure 
9 illustrates this operation. 

The midpoints ' thus corrected are then* summated 
■ during a. second - operation in step*. 3b; , th^ . arrow ^labelled 
. . . 5 -STACK in Figure. 9, illustrating this operation.* %m^ft y hf 
* These two MIG and, STACK operations are specific:, to 
.the second embodiment of, the .invention,. (PSTM migration) . 
- < v T he . following J :,, operations;- .are performed in^^the 
\ m .context of fc each, .of . the j,two embodiments discussed ..(the NMO 
1°, : . 1} oblic^eness >; correction,. as :. ; ^we s ,, have .. .. described/; 
corresponding to the special , case : of a PSTM migration 
with zero migration aperture). ' v :- 
' ; " semblance is : then . calculated for each picking 

time (second operation in step 3a for the NMO correction 
15 only, third operation in step 3b for the PSTM migration) , 
the arrow labelled "semblance" in Figure 9. illustrating 
this opera.tion/v;/ ; / < ' c-\ :k%:-: V./.'.'.::/, T;V 

.operations are implemented for each node, (dtn, i 0 ) • 4 
2 0 , The "automatic bispectral picking" mentioned in 

Figure: 9 corresponds , to picking of maximum semblance * 
1 parameters Jdtn/ x 0 ) for each picking time t 0 (operations ; . ; ' r 
4a,;: 4b in Figures; 8a, 8b respectively) . The < arrow 
labelled MAX illustrates ,the ^ search for;: the, maximum,, 
: 25 semblance. * t ' * r • 

v v Finally, the'pickings of the picked parameters^ (dth, 



'*. Mr-;.-. • V\ -.A. -v. -i-*. 



• ; t 0 ) '. are converted into velocity'"" V ■•'* and':/ ; aneilipt ic ' • 

laws (operations 6a, 6b in Figufes 8&, 8b and arrow 



labelled CONV in Figure 9) . 
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; ^ CLAIMS 

. . , -1. Method . of -determining the, v ^velocity V and 
anellipticity j] -parameters for processing seismic n traces 
in a common midpoint. "(CMP) ; gather including ah^'ahelliptic 
NMO .correction, ^ comprising : J ■ *- v ~ v ^ ^ - T 

5 . > a r preliminary ; step to define a plurality of nodes 

(dtn; '* t 0 ) , the said nbdes^ being ihdi cat i ve of parameters 
dtn and x 0 ^ representing the for " the 

maximum * offset" " and the ' zero "of f set travel time, in 
hyperbolic *" coordinates^ the " said preliminary "step being 
10 followed by *■ > 

- for each node" (dtn, t 0 ) defined in the preliminary 
step, the following steps: 

- for static NMO correction of traces in the CMP 
gather as a function - ! of the values of the said parameters 
15 dtn and x 0 at the node considered, and ... 

for calculating . s ; the , ; semblance . function 
associated " with the said NMO correction for- the node 



considered; * and" - > \ V ' ' " ' ' v 



" for Ve*ach picked"- time t 0 , a ! step including 
mi: 

To ( t 0 ) ) V 



20 determination of . the? maximum semblance node (dtn(t 0 T, 



- and a final step'; to convert the, dtn (t Q ) and, x 0 (to) 
parameters so . as' ;^ to'. , obtain - the^, w velocity ^u^St*)' an 4 ' 
anellepticity T| (t 0 ) laws . . ; ; " ; : '>■'./ " : 

25 2. Method according to claim . 1/ wherein- the nodes 

are defined during the preliminary step in an analysis 
volume (dtn, x 0/ t 0 ) determined by minimum and maximum 
values respectively [dtrimin, dtn^] [t 0m in/ ^omaxl and [t 0 min/ 
tomax] of the dtn 7 x 0 and t 0 parameters . 
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3. Method according to claim 2, wherein, during the 
preliminary step, a corridor [dtn^fto) , dtrwOto) ] , 
[tomin(to), Tomax(to)] for changing dtn and x 0 parameters is 
delimited inside the analysis volume as a" function of 

5 plausible velocity V and anellipticity T| values, / the 
.nodes (dtn, x 0 ) defined for applying 'the" NMO correction 
being then located along the corridor thus ^elimited. 

4. Method according to any one of the preceding 
claims, ' further comprising, for each node. (dtn, \ x 0 ) , a 

0 stacking step of /the corrected seismic traces , following 
the semblance function calculation step. 

5. Method according to claim 4, wherein the stacking 
of * corrected traces is done using only near offset 
traces . 

5 6. Method according to any one of. claims 4 or 5, 

further comprising for each picked time/ and. following 
. the step for determining the maximum semblance node, a 
/step of ' checking that values dtn and x 0 of the maximum/ 
semblance node 5 correspond; to a stacking extreme valued for 

. -7 . 7- . • '• • : Y •' '"' •\'* v .'>. 7 Y • . '• . : ' 7.''.7.';7 ••. ; -,i"> "•••-'>;•..•." -. • Y.t- • ' .;' v -;:' ;;: •••• • 'v. •'-''•"•"■•">'«'.."'" '1' ' ■•' 

'•. • " '• ■ 7 • -7. .' Vv'- >.•'•.• 7v ; -'A •*.-: -"' 77- ; • .' ^ •-7:-7. . ... 7 ■/■:':.■-': 7 : 7' . v ; 7w7*7'H . r7- ." :.:7.7" . 

0 - the same values dtn and x 0 . * \ • < 

■. , 7 . Method r according •,. to any one of the' preceding v 
claims, further comprising a step of selecting and 
adjusting the pickings * obtained, following the step 
implemented for determining; the maximum semblance node 
5 (dtn (t 0 ) , To (t 0 )) for each picked time t 0 , before . the 
conversion step. 

8. Method " according to " claim 7, Wherein /th6^ said 
step of selecting and adjusting the pickings comprises a 
step of only retaining pickings dtn and x 0 for which time 
0 to the highest semblance pickings is greater than a 
predefined value. 
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9. Method according to claim 8, wherein the said 
step of selecting and adjusting the pickings also 
comprises a step, for adjusting the retained pickings dtn 
and t 0 by parabolic , interpolations using values .. about the 

5 ..said picked values - . - %r . t . : •; 

10. Method according to claim 9., wherein the said 
step of selecting and adjustingt pickings also comprises 
a step of eliminating retained and adjusted^ pickings dtn 
and x 0 when it is impossible to calculate the Dix interval 

10 velocities between, the "picking considered: and /higher 
semblance pickings . . , : , • > 

11. Method according to any one of claims 1 to 10, 
wherein the processing applied to seismic, traces is ; . an 
NMO correction process implementing a static correction 

15 CORRnmo- 

12. Method according to claim 11, wherein, during 
the preliminary step, the . NMO corrections , CORRnmo • are 
calculated for all nodes (dtn, To) included in the 
analysis volume and all offsets of processed - seismic 

2 0 traces J ; ' .\ ' }' \ ■/ " 

13 . Method according v to claim 12 , wherein the NMO 

- • ^':Av V^, - "i J^: : '"\ : ' 

correction carried out for each node (dtn, x 0 ) consists of 

applying NMO corrections CORR NM0 calculated during the 
preliminary step. 
25 14;; Method according to any one of claims 11 to 13, 

wherein for a given (dtn, t 0 ) pair, the static NMO 
correction CORRnmo* of a seismic trace with - offset, x is 
carried out according to the following equation: 

CORRnmo (x) « t 0 + \x\ + dtD ^ dtD 1 ^ T °\ 2 in which X max 

V X max 

3 0 represents the maximum Offset in the CMP gather. 
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15. Method . according to any one of claims 1 to 10, 
wherein the processing applied" to seismic traces is' a 
PSTM migration using a static NMO correction CORRp STM . 

16 . Method according - to . claim 15, -- wherein; during 
5 the preliminary step, - .' the NMO., corrections;' CORRp ST M are . 
calculated for all nodes.; {dtn * and t 0 Y' included in the 
analysis volume anci . all migration, .offsets^ inside., the 
y* migration ^ aperture' *'T"- ; ^ r ' ^ 1 * * ^ , *' ' \ '.y 

17. Method according to claim 16, Wherein <* the"" NMO 
10 Vl - correction step ? carried out for each node (dtn : ~and To) v ' 
comprises, for each offset class, application of the said 
NMO corrections CORR PSXM , . calculated .. during . the./ 
- preliminary step on>. all. - midpoints inside the : migration 
aperture. -<\ • : 

15 .18. Method according to claim 17, wherein the NMO 

correction step carried out for;" each npde^ (dtn and t 0 ) 
comprises, for each offset class, the stack of the 
corrected midpoints following application of the said NMO . 
- corrections - CORR PST m. ^ - ^ s 1 v * ''V ? 

20 y 19'/, Method, according, to any : one ■ of /claims ; 15 to 18,' / 

wherein, for a given pair, (dtn arid-- x 0 ) , the 1 ; static "NMO • • \; 
correction CORR PSTM isj carried out according/ to the 
following' equation: ' ' : ' 

Corr PS ™ (x) = -x 0 + jr.*-, dinym + 2Tj(x-x.+l>r t+ f*T*. : <to<M + 2r 0 )(±-x-hY z. 

25< '•>•>■"•• ,/ whete: ' — - -^^y^^-y'^ ;«<<-vvV * - ;.. ';. 

- Xm represents- 'the;: ppordinat:es'''of ' the ''midpoint is ^ 

- x - x ro represents the migration- aperture PSTM, 

- h is the ha.lf source - receiver offset, 

- xmax is the maximum offset and aperture of the 
3 0 migration. 
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" 32 . • 

20. Method according to any one of claims 14 or 19, 
wherein, during the final conversion step, the parameters 
dtn (to) and (x0) are converted; to the velocity law V(tO) 
according to the following equation:.. 



|dm(^ + 2r 0 )^ 



21. Method according to one of t claims 14 , : or 19, 
wherein, - during the final conversion- step> the parameter 
TO (t0) . . is converted tp^ the anellepticity t r\ (to.) ,v law 



according" to r\ = «■ 



^~ -1 



to 



10 : 22 . Method according' 7 to ' claims = 20 and 21; wherein 

parameter dtn is defined with respect to the velocity V 
and anellepticity i\ . according to' the following equation : 



dtn = 




1 + 8T1- ll + 8Tj^ (l + 8r|)V 2 



15 23. Method according; to claim 21, /wherein parameter 

. x0 is defined; according , to anellepticity; ty ;accopiiTig to . 
. the following equation; / - *i - > ; ^ .1 



T 0 = 



l + 8r| 



, 24. Method ' of . characterising a velocity, field for,: . 

2 0 processing seismic data using a gather of seismic . traces ■ 
at common midpoint ■;; : Vherein'V;:- for each. tra.yel tim^ to, for 
a zero offset, a set of parameters dtn and t0 is defined, 
representing the NMO correction for maximum offset, and 
the zero offset travel time respectively, in hyperbolic 

2 5 coordinates . 
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• Determine limi ts of the analysis volume (din, to, to) 

• Calculate NMO corrections CORRmdo for all offsets and 

for all nodes (dtn, t 0 ) 
Calculate corridor limits of the analyses made 




3a 




For each to along the corridor: y i: ; yv, ; 

• v c Application of NMO corrections CORRnmo 

• \ Calcoilate the semblance \/: • \ \ . *. ' • * , 
v Calculate the summation in near offsets- > 



Search for the ma ximum semblance in the corridor and 
the corresponding node (dtn. To) 




Create the dtn(to), To(t6) and semblance (to) series $ 



Increasing sort of the semblance series (to) . 

Reject pickings too close to picldngsl^ with strong 

semblance # ' ;<* y . V 

Adjust selected values (dtn, to) by parabolic interpolation 

Reject pairs (dtn, to) for which the Dix interval velocities 

with pairs (dtn, to) with the strongest semblance, are 

unacceptable 
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FIG. 8 b 
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lb 



Determine limits of the analysis volume (dta, to, to) 
Calculate NMO corrections CORRpstm for all nodes (dtn, 
T 0 ) and offsets inside the migration* aperture 
Calculate corridor limits of the analyses made 





For x each offset class: 0 ; . ; t - 4 { v \ \ . x 

• ' Application of NMO corrections CORRpstm onfall midpoints 
inside the migration aperture, along the corridor ' 

• Summation of corrected midpoints along the comdor >; - 
For each ^ along the .corridor: 0 V 1 , 

• Calculate the semblmce^ : . V., , . v \. v - ■ 

• Calculate the summation in offsets ,. 



Search for the maximum semblance in the corridor and 
the corresponding node (dtn, T 0 ) 




Create the dtn(to)f T 0 (to) and semblance(to) series 




v-. Increasing sort of the semblance series (fy) ^ 

• Reject pickings too close to pickings with strong 
: : v;; semblance > /<;■; 1 1 * 

• Adjust selected values (dtn, t 0 ) by parabolic interp olation 

• Reject pairs (dtnf x 0 ) for which the Dix interval (dtn, t 0 ) 
with the strongest semblance, are unacceptable 
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